Optimal conditions for exploring high-density baryonic matter 
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Using simple parametrizations of the thermodynamic freeze-out parameters extracted from the 
data over a wide beam-energy range, we reexpress the hadronic freeze-out line in terms of the 
underlying dynamical quantities, the net baryon density pB and the energy density e, which are 
subject to local conservation laws. This analysis makes it apparent that ps exhibits a maximum as 
the collision energy is decreased. This maximum freeze-out density has p — 400 — 500 MeV, which 
is above the critical value, and it is reached for a fixed-target bombarding energy of 20 — 30 GeV/A. 



PACS numbers: 25.75.-q, 05.70.-a, 64.70.-p, 64.90. +b 

In recent years, it has become abundantly clear that 
the collision energy plays an important role in determin- 
ing the properties of the final state in relativistic heavy 
ion collisions. While the extracted freeze-out tempera- 
ture generally increases monotonically with the collision 
energy, the corresponding net baryon density exhibits a 
more complicated behavior: It initially increases steadily 
as the collision energy is raised but ultimately, at suffi- 
ciently high energies, it decreases as a result of the nu- 
clear transparency. Therefore, there is a certain beam 
energy for which the resulting baryon density at freeze- 
out acquires a maximum value. In this paper, we discuss 
this optimal beam energy on the basis of the most up-to- 
date results on the properties of the hadronic freezeout. 

For this purpose, we employ statistical considera- 
tions which have been quite successful in accounting for 
the observed yields of most hadronic species (see Refs. 
B II II II II . for example). In fact, the evidence has 
become truly overwhelming in recent years that chem- 
ical equilibrium is established in relativistic heavy ion 
collisions throughout the entire range of collision ener- 
gies . The resulting picture shows that the freeze-out 
temperature T increases steadily with the collision en- 
ergy as the corresponding baryon chemical potential ps 
decreases steadily towards zero. 

However, it is not always the most convenient to work 
with the thermodynamic variables T and (is- This is 
particularly true within a dynamical context, since they 
are not subject to any conservation laws, in contrast to 
the more basic dynamical variables, namely the energy 
density e and the net baryon density ps- Furthermore, 
when a first-order phase transition is present, T and ps 
become non-monotonic functions of e and pb throughout 
the associated phase-coexistence region of the phase dia- 
gram (as does the pressure p) . This feature in turn causes 
the inverse functions to be multivalued, a clearly incon- 
venient situation. It is therefore of interest to reexpress 
the thermodynamic freeze-out information in terms of 
the dynamical quantities. The resulting representation 
also serves to more clearly bring out the fact that the 
freeze-out density exhibits a maximum value, a feature 
that may be important in the planning of experiments 
that seek to explore compressed baryonic matter. 



In the grand-canonical ensemble of hadron species i, 
the partition function factorizes into separate contribu- 
tions for each specie, Z = HiZi, with 
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where plus is for bosons and minus is for fermions. The 
hadron mass is nii and g-i=2Ji+l is the spin degeneracy. 
The fugacity associated with the hadron specie i is 

Ai(T, {p}) = c (kbB,+h q q 1+ hsS,)/t ^ 

where Bi, Qi, and Si denote its baryon number, electric 
charge, and strangeness, respectively. The ensemble is 
characterized by the temperature T and the three chem- 
ical potentials {p,} = {p,B, I^Q, Ms}- Finally, V denotes 
the enclosing volume of the system (which plays no role 
in the analysis since only yield ratios are considered). 

The corresponding spatial number density of a given 
specie is then readily obtained, 
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In the first line, the factors are arranged such that the 
terms in the sum are regular in the small-mass limit, 
mi — > 0, where K 2 diverges, and the second line ex- 
hibits the leading term representing the classical result. 
Since the sign of the second term in Q depends on the 
quantum-statistical nature of the specie, the leading term 
overestimates fermion densities while it underestimates 
boson densities. For all the parameter values employed 
in the yield estimates, the first term in the expansion, Eq. 
(0J, is an excellent approximation for all baryon species, 
being accurate to within a few per mille. For kaons it is 
off by a few per cent, and even the pion densities are un- 
derestimated by at most 10%. The first term |(2J would 
thus be a quite reasonable approximation for rough esti- 
mates of the yields. 

Once we know the number density of each specie, it 
is straightforward to calculate the corresponding baryon, 
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charge, and strangeness densities. In particular, the net 
baryon density is 

p B (T,{p}) = = 5>^(T,M) . (5) 



Furthermore, the total energy density is e = J^. £j, where 
the contribution by a particular specie i is 
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where the last line is the classical result. 

As it turns out, it is possible to fit the observed yield 
ratios reasonably well with the statistical model, by ad- 
justing the Lagrange multipliers at each particular colli- 
sion energy 0,H,US0- I n this analysis, T and ps are 
treated as free parameters, while the values of ps and 
Pq are fixed by requiring overall strangenss neutrality, 
(S) = and that the ratio of the net charge to the net 
baryon number be equal to that of the collision system, 
e.g. (Q) = 0.4(B) for Au+Au. As discussed in Ref. @, 
the effect of the small change in this ratio in going to 
Pb+Pb is negligible. The extracted values of the freeze- 
out temperature can be approximately represented as , 
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where T and pb are expressed in MeV. 

Using the above relationship JSJ), we have calculated 
the corresponding freeze-out values of the net baryon 
density ps (Eq. 10) and the energy density e (Eq. |J§}). 
The result is displayed in Fig. ^ Before turning to the 
calculated results, we note that the lower-right section of 
the pb — £ plane is not physically accessible, since a given 
net baryon density ps gives rise to a certain minimum 
energy density. (The corresponding curve, e min (pB), i s 
simply the pressure at zero temperature and is therefore 
occasionally, and somewhat misleadingly, referred to as 
the equation of state.) In the ideal-gas scenario, where 
binding and compression effects are absent, this lower 
bound is given by e m in = ttinPb, which is indicated on 
the figure. 

In order to illustrate the importance of adjusting pq 
and ps to ensure average conservation of charge and 
strangeness (see above), we have also calculated those 
values of pb and e that would result if pq and ps were 
taken to vanish. As can be seen from Fig.^ the failure to 
adjust pq and ps increases the values of both densities 
by amounts that are epsecially significant in the region 
of maximal density. These results were calculated both 
clasically (first terms only) and quantally (all terms) in 
order to illustrate how unimportant this distinction is for 
the present analysis. 
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FIG. 1: The hadronic freeze-out line in the pB ~ e phase 
plane as obtained in the statistical model with the values of 
Pb and T that have been extracted from the experimental 
data in Ref. 0- The curves on the right have been calculated 
for pq = us = using either quantal (solid) or classical 
(dashed) statistics, while the curve on the left employs values 
of pq and ps that have been adjusted to ensure {S) = and 
(Q) = 0.4(B) for each value of pb- Also indicated are the 
beam energies (in GeV/A) for which the particular freeze-out 
conditions are expected at either RHIC (total energy in each 
beam), starting at 100+100 and going down to 2+2, or FAIR 
(kinetic energy of the beam for a stationary target), starting 
at 5 and going up to 40. as based on fits to existing data |g. 
The triangular area corresponds to energy densities below the 
minimum required at the given net baryon density, s = mjvps 
(ignoring binding and compression), and is thus inaccessible. 



At the highest energies, freeze-out occurs for a negli- 
gible value of the chemical potential ps , (and hence the 
net baryon density ps is practically zero) and the energy 
density e is nearly one half GeV/fm 3 . As the collision 
energy is lowered, p increases rapidly while T initially re- 
mains fairly constant but gradually begins to drop. Then, 
in the range of T = 140 — 130 MeV, the freeze-out line 
{Pb,s) exhibits a backbend and approaches the origin. 
The resulting maximum value of the net baryon density 
at freeze-out is about three quarters of the normal nu- 
clear saturation density of po ~ 0.15 fm -3 . 

The value of the baryon chemical potential ps ex- 
tracted from the data decreases monotonically with the 
collision energy and can be parametrized on a simple 
form 
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where ps as well as the NN CM energy y/s are expressed 
in MeV. By use of this result, it is possible to attach 
beam energies along the freeze-out curves shown in Fig. 
^ This has been done for both collider experiments, such 
as those being carried out at RHIC, or for a fixed target, 
as has been done at the SPS and the AGS and is being 
planned at FAIR. 

We note that if this region of maximum freeze-out 
density were to be explored with RHIC, the appropri- 
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ate beam kinetic energies would be about 2 — 4GeV/A, 
corresponding to s = (6 — 10GeV/^4) 2 , which would be 
a rather challenging proposition. By contrast, it ap- 
pears that the region would be well within reach of FAIR 
which is planned to bombard fixed targets with heavy-ion 
beams having kinetic energies of 5 — 40 GeV/A. 

Finally, according to the most refined lattice QCD 
studies H 0, 0, 0, 0| , it appears that the region of 
highest freeze-out density lies well beyond the critical re- 
gion where the transformation between the quark-gluon 
plasma and the hadronic resonance gas changes from be- 
ing merely a crossover to being a genuine first-order phase 
transition. This fact is helpful to the prospect of probing 
this phase transition by means of heavy-ion collisions. 

Let us briefly summarize: Analyses of experimentally 
obtained hadronic yield ratios at a variety of collision 
energies have shown that the data can be well repro- 
duced within the conceptually simple statistical model 
that describes an ideal hadron resonance gas in statis- 
tical equilibrium. Furthermore, the extracted freeze-out 
values of the temperature and the baryon chemical po- 
tential exhibit a smooth and monotonic dependence on 



the collision energy and can be simply parametrized. We 
have used these results to examine how the freeze-out 
appears when represented in terms of the basic baryon 
and energy densities, rather than chemical potential and 
temperature. These quantities are more basic and they 
are of more direct relevance to the collision dynamics, 
because they are subject to corresponding conservation 
laws. 

We have found that the freeze-out value of the net 
baryon density exhibits a maximum as the collision en- 
ergy is being scanned, thus suggesting that there may be 
an optimal collision energy (range) for the exploration of 
compressed baryonic matter. In a fixed-target configura- 
tion, this optimal beam kinetic energy is 20 — 30 GeV/A, 
which appears to be ideal for the planned FAIR at GSI, 
and it may be accessible in the low-energy campain now 
being planned at RHIC. 
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